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Abstract 

It is shown by numerical calculations that the convoluted forward pomeron propagator in the 
external field created by a solution of the Balitski-Kovchegov equation in the nuclear matter vanishes 
at high rapidities. This may open a possibility to apply the perturbative approach for the calculation 
of pomeron loops. 

1 Introduction 

In the QCD, in the limit of large number of colours, strong interaction at high energies is mediated 
by the exchange of BFKL pomerons, which interact via their splitting and fusion. In the quasi- 
classical approximation for photon (hadron)-nucleus scattering the relevant tree (fan) diagrams are 
summed by the well-known Balitski-Kovchegov (BK) evolution equation [1, 2, 3]. For nucleus- 
nucleus scattering appropriate quasi-classical equations were derived in [4, 5]. In both cases 
pomeron loops were neglected. This approximation can be justified if the parameter 7 = A exp Ay 
is small, with y the rapidity and A and A the pomeron intercept and triple pomeron coupling. 
Then for a large nuclear target, such that A 1 ^ 3 ^ ~ 1, the tree diagrams indeed give the dominant 
contribution and loops can be dropped. However with the growth of y the loop contribution becomes 
not small and this approximation breaks down. 

Direct calculation of the loop contribution seems to be a formidable task for the non-local 
BFKL pomeron. Simplest loops have been studied in several papers for purely hadronic scattering 
[6, 7, 8]. In particular in [8] it has been found that pomeron loops become essential already at 
rapidities of the order 10-i-15. They shift the position of the pomeron pole to the complex plane and 
thus lead to oscillations in cross-sections. However with the growth of energy loop contributions 
begin to dominate and one needs to sum all of them. There have been many attempts to do this 
in the framework of the so-called reaction-diffusion formulation of the QCD dynmaics and the 
following correspondence with the statistical approach [9, 10, 11, 12, 13, 14] (see also a review 
[15] and references therein). Unfortunately concrete results could be obtained only with very 
crude approximations for the basic BFKL interaction and the stochastical noise in the statistical 
formulation. The conclusions of different groups are incomplete and contradictory. So in [14] it was 
found that the geometric scaling following from the BK equation was preserved with loops taken 
into account, although going to the black disc limit was much slower. On the contrary in papers 
based on the analogy with statistical phyiscs (see [13, 15]) it was argued that the BK scaling was 
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changed to the so called diffusive scaling (with an extra y/y in the denominator of the argument) 
but the speed of achieving the black disk limit was essentially unchanged. 

In our previous study of pomeron loops [16] we considered a much simpler model with the 
local supercritical pomeron in the Regge-Gribov formalism. Instead of trying to solve the model 
for the purely hadronic scattering we considered the hadron-nucleus scattering and propagation 
of the pomeron inside the heavy nucleus target. Moreover to avoid using numerical solution of 
the tree diagrams contribution with diffusion in the impact parameter, we concentrated on the 
case of a constant nuclear density which allowed to start with the known analytical solutions. We 
have found that the nuclear surrounding transforms the pomeron from the supercritical one with 
intercept e > to a subcritical one with the intercept — e. Then Regge cuts, corresponding to loop 
diagrams, start at branch points located to the left of the pomeron pole and their contribution is 
subdominant at high energies. As a result the theory aquired the properties similar to the Regge- 
Gribov with a subcritical pomeron and allows for application of the perturbation theory. In [16] we 
expressed our hopes that a similar phenomenon might occur in the QCD with BFKL pomerons. 

In this note we demonstrate that such hopes are possibly founded. We consider the pomeron 
propagator in the external pomeron field created inside the nucleus and give arguments that, similar 
to the local Regge-Gribov case, it vanishes at large rapidity distances. We stress that at present we 
are unable to give the full proof for this behaviour. Our study is based on numerical calculations. 
This makes us to choose a relatively small subset of initial conditions out of the complete set 
necessary for the study of the pomeron propagator. Moreover, due to technical difficulties, in 
this note we restricy ourselves to the much simpler pomeron propagator in the forward direction. 
Our numerical results show that, with the chosen set of initial conditions, this forward propagator 
vanishes at large rapidity distances. This result is insufficient for the study of loops, where non- 
forward propagator are involved. However it can be applied for double inclusive cross-sections in 
the nucleus- nucleus scattering, in which only forward propagators are important. 

2 Main equation 

We are going to study the behaviour of the BFKL pomeron propagator in the external field, 
generated by a solution of the BK equation. We consider the simplified case of the nuclear matter, 
when the dependence on the impact parameter b is absent. This propagator corresponds to a sum 
of diagrams shown in Fig. 1. 

As mentioned, in this note we restrict ourself to the forward propagator, P(y, x, x') when the 
transferred momentum is zero. It depends on rapidity y and two 2-dimensional coordinate vectors x 
and x' corresponding to the initial and final distances between the reggeized gluons in the pomeron. 
At y = we have P(y, x, x') = X7~ 4 5 2 (x — x'). Since the study is only possible numerically, to avoid 
using this singular initial condition, we shall consider a convolution of P(y, x, x') with an arbitrary 
initial function V 4 ip(x) 



This convolution satisfies the same equation as the propagator itself but at y = we have 



Obviously properties of the propagator can be studied taking a full set of functions ip(x). 




(1) 



P(y = 0,x) =1>(x). 



(2) 
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Figure 1: Diagrams summed into the pomeron prpagator in the nuclear matter 



The equation for P(y, x) can be conveniently obtained from the BK equation for the sum of fan 
diagrams. In the forward direction this sum Q(y,x) satisfies a non-linear eqiation 



d$(y,x) 



OL f x^ / \ 

— / d 2 x l ^^[^(y,x 1 ) + &(y,x 2 ) - &(y,x) - $(y,x 1 )Q(y,x 2 )), (3) 

Z7T J X\X2 V 1 



ot s N c 



dy 

where standardly 

a = — 

The equation for the convoluted propagator P(y, x) in the presence of nuclear medium is obtained 
when one of the $ in the non-linear term in (3) is substituted by a particular solution of (3) with 
a given boundary condition. The equation thus obtained is 

dP(y,x) _ a f 2 x 



dy 



I d 2 x 1 ^(p(y,x 1 ) + P(y,x 2 ) - P{y,x) -2$(y,x 1 )P(y,x 2 j). (4) 

Z7T J X\X^ \ / 



Note that the initial condition for <&(?/ = 0,x) = &o(x) is fixed by the properties of the nuclear 
medium, whereas the initial condition (2) for P(y, x) is arbitrary, since we we are interested in the 
propagator in a given nuclear surrounding. 

Equation (4) is a linear equation for P(y, x) in contrast to the BK equation. At y — > oo 
<&(y,x) —7- 1 independent of the chosen initial condition. One may think that at y — > oo the 
behaviour of P(y, x) can be derived from the asymptotic equation 



dP(y,x) 



dy 



y— >oo 



I d 2 x 1 ^(p{y,x 1 )+P(y,x 2 )-P(y,x)-2P(y,x 2 )) 



a / , 2 x 
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= -P(y,x)— / d xi-5-2. (5) 
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However the integral on the right-hand side has become divergent (although it converges at finite 
y). This means that the limit y — > oo is more delicate and cannot be taken under the sign of integral 
over x\. And indeed we shall see by numerical calculation that the behaviour of the solution at 
y — > oo is not solely determined by the limiting value of $(y,x) but depends on its behaviour at 
finite y. 

For numerical studies both the BK equation and linear equation (4) in the momentum space 
are more convenient. Introducing 

(f>(y,x) = 2—, p{y,x) = — (6) 

x z x z 

and then passing to the momentum space we obtain the following equations for <p(y, k) and p(y, k) 

d(f)(y, k) 



and 



dy 



dp(y, k) 



= -a(H BFKL 0(y, k) + 4? (y, kfj (7) 
-o.{h B fkl + 2<t>(y, k)^p(y, k), (8) 



dy 
where 

Hbfkl = In k 2 + In x 2 - 2(^(1) + In 2). (9) 

To study the behaviour of the propagator in the external field <p one has to solve this pair of 
equations with the initial conditions 

4>(y, k) y=0 = 4> (k), p(y, k) y=0 = p (k), (10) 

with some fixed cj>o and for a complete set of function po(k). 



3 Numerical studies 

We have set up a program which simultaneously solves the pair of equations (7) and (8) for a given 
pair of initial conditions (10). For the BK evolution we have fixed the initial condition as 

cf) (k) = -(l/2)Ei(-fc 2 /0.3657) (11) 

used in our previous calculations. The behaviour of 4>(k) with k 2 at different values of the scaled 
rapidity Y = ay = 2, 4, 6, 8 and 10 is shown in Fig. 2. (Note that the maximal value of the scaled 
rapidity Y = 10 corresponds to the natural rapidity of order 50). 

For the BFKL evolution in the external field (f>, in the first run (A), we have taken the same 
form of the initial condition but with a variable slope 

po(fc) = -(l/2)Ei(-A; 2 /a). (12) 

We have performed calculations for a = 0.2, 0.6, 1.0, 1.4 and 1.8. In the second run (B) the initial 
condition was taken with extra powers of k 2 . 

p (k) = -(l/2)A; 2n Ei(-A; 2 /0.3657) (13) 

with n = 0,1,2,3 and 4. 
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Figure 2: Solution (f)(y, q) of Eq. (7) as a function of q 2 for different of Y = ay. Curves from 
bottom to top correspond to Y = 0, 2, 4, 6, 8 and 10. 

In all cases the behaviour of the solution p(y, k) was found to be universal. At large enough y 
the solution becomes independent of k 2 up to a certain maximal k'^ nax (y), starting from which it 
goes to zero. Roughly 

p{y,k)~A{y)9{k 2 max -k 2 ). (14) 

As y grows A(y) goes to zero and km ax (y) goes to infinity. . So on the whole the solution vanishes 
as y — > oo, its x dependence tending to 5 2 (x). 

We illustrate this behaviour in Figs. 3 and 4, in which we show the solution p(y, k) for run A 
with a = 1.0 and run B with n = 2 as a function of k 2 . One observes that although the values of 
p(y, k) for the two cases are different, their behavior with y is the same: they vanish as y — > oo. 

This unversality is especially obvious if one calculates the slope A(y, k) of the y-dependence of 
p(y, k) at fixed k presenting 

p(y,k) oc e YA(y ' k \ (15) 

It turns out that at Y > 1 the slope A(y, k) is independent of k and identical for all considered 
cases (run A with all studied a and run B with all studied n). Its smooth behaviour with Y is 
shown in Fig. 5. One observes that starting from Y = 5 the slope becomes negative indicating that 
the solution goes to zero at Y » 1. 

One has to take into account that in the external field <f)(y, k) depending on rapidity the pomeron 
prpagator ceases to depend only on the rapidity difference. Rather the initial and final rapidities 
become two independent variables. To see what influence it has on the behaviour of the propagator 
at large rapidities we varied the initial rapidity y = yo for the evolution of p(y, k), leaving unchanged 
the initial rapidity y = for the evolution of 4>(y, k), which is the rapidity of the nucleus. One finds 
that although at initial stages of evolution the behavior of p(y, k) strongly depends on the value 
of yo, at higher rapidities this behaviour is essentially the same for any yo, namely the convoluted 
propagator goes down with rapidity with the slope independent of yo- This is illustrated in Figs. 
6 and 7 which show results for yo = 3/a (Yq = 3) In Fig. 6 we show the solution p(y,k) for run 
A with a = 1. One observes, that although absolute values of p(y, k) are quite different from the 
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Figure 3: Solutions p(y, q) of Eq. (8) as a function of q 2 for different Y = ay for run A with a = 1. 
Curves which start to fall at higher g 2 correspond to higher Y = 0, 2, 4, 6, 8 and 10 
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Figure 4: Solutions p(y, q) of Eq. (8) as a function of g 2 for different Y = ay for run B with n = 2. 
Curves which start to fall at higher q 2 correspond to higher Y = 0, 2, 4, 6, 8 and 10 
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Figure 5: Slope A of the y-dependence of the solutions p(y,q) of Eq. (8). The slope is the same 
for all initial conditions (12) and (13). 



case yo = shown in Fig 3, the behaviour with the growth of rapidity is the same. It is especially- 
clear from the values for the slope A shown in Fig. 7 together with those for the case y = (Fig. 
5). Again at the initial stage of the evolution the behavior with Y$ = 3 is quite different from that 
with Yq = 0. However at higher rapidities the values for the slope become the same. 

It is remarkable that this behaviour takes place only with cj> given by the exact solution of the 
BK equation. Taking an approximate form 

<S>{y,x)~l-e- Q2 ^ x2 , (16) 

where the "saturation momentum" Q 2 (y) ~ e 2may we obtain an equation for p(y,k) in the mo- 
mentum space 

^ H-nm)] <»> 

Taking for simplicity the initial condition po(k) = 4>o(k) at y = we get the solution shown in Fig. 
8. One observes that at large y the solution acquires the same form (14) where however A(y) grows 
with y: 

A(y) ~ e 1AY (18) 

This implies that with the approximate form (16) of <fr the final solution P(y, x) in the coordinate 
space behaves in a singular manner at y — > oo. Effectively 

P(y,x) y ^ 00 ^e 1A y X 2 5 2 (x) (19) 

and it is impossible to say that it vanishes in this limit. 



4 Conclusions 

We have studied numerically the BFKL pomeron forward propagator in the external field created 
by the solution of the BK equation in the nuclear matter. We have found that for more or less arbi- 
trary set of initial conditions the convoluted propagator vanishes at large rapidities, its coordinate 
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Figure 6: Solutions p(y, q) of Eq. (8) as a function of q 2 for different Y = ay for the initial value 
taken at Y = 3 according to run A with a = 1. Curves which start to fall at higher q 2 correspond 
to higher Y = 4, 6, 8 and 10 



Figure 7: Slope A of the ^-dependence of the solutions p(y, q) of Eq. (8) with the initial condition 
put at Y = (upper curve) and at Y = 3 (lower curve) according to run A with a = 1. 
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Figure 8: Solutions of Eq. (17) as a function of k 2 for different y. Curves from bottom to top 
correspond to Y = ay = 0, 2, 4, 6, 8 and 10. 

dependence tending to the <5-function. This gives reasons to believe that the forward propagator 
itsef vanishes at large rapidities in the nuclear background. This result follows only with the field 
being the exact solution of the BK equation. 

Our results are obviously insufficient for the calculation of pomeron loops, which requires the 
non-forward pomeron propagator. However they can be directly applied to the study of double in- 
clusive cross-section for gluon jet production in nucleus-nucleus collisions. This highly complicated 
problem is left for future investigation. 
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